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ABSTRACT 

We study the leading discretization errors for staggered fermions by first constructing the con- 
tinuum effective Lagrangian including terms of 0{a?), and then constructing the corresponding 
effective chiral Lagrangian. The terms of O(a^) in the continuum effective Lagrangian completely 
break the S'C/(4) flavor symmetry down to the discrete subgroup respected by the lattice theory. 
We find, however, that the O(a^) terms in the potential of the chiral Lagrangian maintain an S'0(4) 
subgroup of SU{4:). It follows that the leading discretization errors in the pion masses are 5*0 (4) 
symmetric, implying three degeneracies within the seven lattice irreducible representations. These 
predictions hold also for perturbatively improved versions of the action. These degeneracies are 
observed, to a surprising degree of accuracy, in existing data. We argue that the 50(4) symmetry 
does not extend to the masses and interactions of other hadrons (vector mesons, baryons, etc), nor 
to higher order in a^. We show how it is possible that, for physical quark masses of O(a^), the new 
SO (4) symmetry can be spontaneously broken, leading to a staggered analogue of the Aoki-phase 
of Wilson fermions. This does not, however, appear to happen for presently studied versions of the 
staggered action. 
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1 Introduction 



Staggered fermions are one of the standard discretizations used in lattice simulations. 
Their major advantage is that they retain a remnant of the continuum chiral symmetry group 
at non-zero lattice spacing. Their major drawback is that the continuum flavor symmetry 
{SU{4:) for a single staggered fermion) is almost completely broken on the lattice — only a 
discrete subgroup remains. Although only an effect of O(a^), flavor breaking turns out to be 
numerically significant at typical lattice spacings. 

In this paper we show that, in the Goldstone boson sector, flavor symmetry is partially 
restored in the chiral limit. The discrete lattice flavor group enlarges to an 5*0(4) subgroup 
of the continuum SU{4). This implies relations between multi-pion interactions at very 
low momenta. The simplest of these concerns the pion masses: seven lattice irreducible 
representations (irreps) collapse into four 5*0(4) irreps in the chiral limit, implying three 
degeneracies. This result applies to both the original form of staggered fermions, and to 
various improvements that have been recently suggested H, |, |, ^. It turns out that these 
degeneracies are present both in data from a few years ago 0, and from recent simulations [Q. 

Our derivation of these degeneracies applies only to the leading order discretization errors. 
We find that they are violated by terms of 0(a'^). We have also studied the extension of 
our result to the masses and interactions of hadrons other than the pions — we find that the 
partial flavor symmetry restoration does not extend to the masses and interactions of the 
vector mesons, and we argue that this is also true for any particle which remains massive in 
the chiral and continuum limits. In other words, for such particles all lattice irreps should be 
split by terms of O(a^), even in the chiral limit. Present numerical results for vector mesons 
and baryons are either of insufficient accuracy or do not consider enough irreps to test this 
prediction. 

In the bulk of this article we explain how the enlarged flavor symmetry is derived. The 
method is modeled on that used to study flavor and parity breaking with Wilson fermions [0. 
While straightforward in principle, the application of the method to staggered fermions is 
considerably more complicated. One begins by determining all the lattice operators which 
are required for O(a^) improvement. This has largely been done by Luo P], but he missed 
some operators, and we complete the task in App. The second step is to write down the 
effective continuum action which describes the lattice theory including errors proportional 
to a^, and determine the symmetries of its various components. This is done in Sec. ^ 
A side benefit of this analysis is seeing explicitly how the lattice symmetry group emerges 
in continuum language. One then writes down the effective chiral Lagrangian describing 
the interactions of the pseudo-Goldstone bosons in the effective continuum theory, again 
including the terms. This is worked out in Sec. § (with details relegated to App. P), 
and it is at this stage that the symmetry restoration emerges. The implications of partial 
symmetry restoration for pion masses are presented in Sec. ^, and are compared to the results 
of numerical simulations. In Sec. |^ we explain why the partial symmetry restoration does 
not extend to vector mesons. 

It turns out that the enlarged flavor symmetry which emerges in the chiral limit at 
O(a^) can also be spontaneously broken in that same limit. Whether this happens depends 
on the sign of various unknown coefficients in the chiral Lagrangian describing the lattice 
theory. This phenomenon is analogous to the breaking of flavor in the Aoki-phase for Wilson 
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fermions |^. The analogy is not precise, however, because the symmetry that is being 
broken here is not an exact symmetry of the lattice theory, but is instead violated by terms of 
0{a'^). Thus there are only pseudo-Goldstone bosons in the broken phase. The new phase, if 
it occurs, is very narrow — in terms of the bare lattice quark masses it has a width of O(a^). It 
appears that this phase is unlikely to occur for those versions of staggered fermions presently 
in use, but it may occur for other improved versions. In Sec. |^ we give a brief discussion of 
the properties of the new phase. 

We close with some comments on the practical implications of our result. 

A preliminary report on this work was presented in Ref. At that time we suggested 
that the enlarged symmetry would apply also to hadrons other than the pions, a possibihty 
we now think is incorrect. 



2 The Continuum Effective Action 

The first step in our analysis is to write down the effective continuum action describing 
the interactions of quarks and gluons with momenta much below the cut-off, p -C A ~ 1/a 



111 . Integrating out modes near the cut-off introduces higher dimension operators with 
coefficients suppressed by explicit powers of a. The constraint on these operators is that 
they must have the same symmetries as the underlying lattice action. We wish to include 
terms suppressed by up to two powers of a, and thus must determine all allowed operators of 
dimension up to and including six. This we do by writing down the allowed lattice operators 
and then matching them onto continuum operators at tree level. The result has the form 

Scs = 84 + 0^3^ (1) 
= J d^x{C4 + a^Ce) . (2) 



The leading term is 



with 



A = ^TrF^uF^, + + m)Q , (3) 



^ = E(7m®1)^m, (4) 

in a (spin ® flavor) notation explained in more detail in App. 0. The physical quark mass 
is proportional to the bare staggered quark mass m(/i)a = ^'^(5'^, ln(a/x))mo(a). There are 
no terms linear in a because there are no operators of dimension five consistent with all 
the lattice symmetries. This is well known for the gauge theory, and has been shown in 
Refs. O, O] for staggered fermions. 



There are three types of contribution to S^. gluonic operators, fermion bilinears and 
four-fermion operators 



The gluonic terms are [|14 



~ Tr {D^F^^DpFp,) + ^ Tr (D^F^.D^F,,) + J] Tr {D^F^.D^F^,) , (6) 
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where the symbol ~ imphes that each of the operators appears with a coefficient which 
depends on and In a. The ffist term in eq. (P) can be absorbed into a redefinition of gluon 
field, which leaves only two independent terms in Cf^°. 

The enumeration of the fermionic operators is described in App. along with the nota- 
tion used to refer to them. Most have been previously listed by Luo ||^. There are 8 fermion 
bilinears 

+ mQ{ipfQ + mQDlQ + m^QjpQ + m'QQ , (7) 

and 24 four-fermion operators. The latter we divide into two classes: 

£f (A) ~ [S ^ A] + [S ^V] + [Ax S] + [V ^ S] + 

+ [P xV] + [P X A] + [V X P] + [Ax P] + 

+ [T X V] + [T X A] + [V X T] + [A X T] + 

+ {[S X 5] - [P X P]} + {[S X P] - [P X S]} + 

+ {[S X T] - [P X T]} + {[T X S] - [T X P]} + 

+ {[V xV]-[Ax A]} + {[V X A]-[AxV]} (8) 

(with compound operators enclosed in curly braces), and 

rf^^) ~ [T, X 1^^] + [T^ X A^] + [K^ X + [A^ X + 

+ {([l^, X V,] - [A, X A,]) - \{[V X V] - [A X A])} 

+ {{K X ^m] - [^M X V,]) - X A] - [A X V])} . (9) 

The essential point of the notation is that the first letter indicates the rotation property of 
the bilinears within the operator, while the second indicates the fiavor. Thus in [S x T] both 
bilinears are scalars and have "tensor" fiavor, i.e. a fiavor matrix of the form ^^^j, = . 

Several of these operators are redundant, i.e. can be absorbed into the others by a field 
transformation. Although we could use this to reduce the number of operators, we choose 
not to do so, since it is little extra work to keep track of all the operators. In this way we do 
not have to worry about the subtleties of determining redundant operators with staggered 
fermions |p. 

As noted in the introduction, a variety of improvements to the staggered fermion action 
are being tested in numerical simulations. All of these use perturbation theory, and thus 
reduce, but do not remove, the errors proportional to a?. This implies that the form of the 
continuum effective action remains the same as for unimproved staggered fermions, although 
the size of the numerical coefficients multiplying the o? terms is presumably reduced. 
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Term in action [Flavor] x Rotation symmetry 



(m = 0) [SU{A)l X St/(4)^J X 50(4) 
^4 {m ^ 0) [SU{A)] X 50(4) 



Qglue 


[5?7(4)i X 5?7(4)fi] X SW4 


^bilin (^^ ^ 


[SU{A)l X 5?7(4)^] X SW^ 


^bilin _^ 


[SU{A)] X 5W^4 


^FF(A) 


[U{1)a X r4 X 50(4)] X 50(4) 


qFF(B) 


U{1)a X (r4 X 5W^4,diag) 



Se{m = 0) t/(l)A X (r4 X 5l^4,diag) 

56(m ^0) T,x 5iy4.diag 

Table 1: The flavor (including cliiral) and rotation symmetries respected by various terms 
in the effective action. In the last three lines flavor and rotation symmetries are intertwined. 
See text for detailed explanation. 

2.1 Symmetries of the effective action 

For the subsequent analysis, it is crucial to understand the symmetries that are maintained 
by the various terms in the effective action. 

The leading order action 54 is invariant under Euclidean translations, rotations and 
reflections, and also under charge conjugation C, fermion number U{l)v, and the flavor 
group SU{4). In the massless limit the flavor symmetry enlarges to the chiral symmetry 
5t/(4)i X SU{4:)r. In general, we expect 56 to break this symmetry down to that of the 
underlying lattice theory, because we have included in 56 all operators consistent with the 
lattice symmetries. While it is possible that some symmetries will not be broken until higher 
order in a, this turns out not to be the case here. The only exception is continuous translation 
invariance, which is a symmetry of Ses but is broken by the lattice to a discrete subgroup. 
This symmetry is effectively restored because we are considering only low-momentum modes. 

We first note that fermion number, charge conjugation and spatial inversion (parity) are 
not broken by 56- This is obvious for fermion number. Charge conjugation does flip the sign 
of vector and tensor bilinears, but this sign cancels since these bilinears always appear in 
pairs. Parity also flips the sign of some bilinears, but again the overall sign cancels. 

This leaves us to consider the breaking of flavor and rotation symmetry. In Table ^ we 
display the symmetries respected by each of the terms in the effective action. The notation is 
as follows: 5W4 is the hypercubic subgroup of the Euclidean rotation group 50(4); r4 is the 
Clifford group with four generators which square to the identity; X indicates a semi-direct 
product; and U{1)a is the axial symmetry 

Q exp (i^A(75 ® 6)) Q, Q^Qexp {iOAil^ ® 6)) , (10) 
which is the continuum equivalent of the lattice axial symmetry. The meaning of 51^4,diag 
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will be described below. 

We now explain the results in the Table, commenting on their significance as we proceed. 
The gluonic action S*!'"*^ does not affect the flavor symmetry, but the third term in eq. (^) 
breaks the 50(4) rotation symmetry down to its hypercubic subgroup 5*^4, because the 
index fi is repeated four times. The same holds true for Sg''™: all terms in eq. (|^ are flavor 
singlets, but the fourth term breaks 5*0(4) — > SW^. 

Flavor symmetry is broken only by the four-fermion terms. This occurs in two stages. 

FF( A.) 

All terms in have their Lorentz and flavor indices contracted separately, e.g. 

[Txv]^Y. E^(^/)(7m^ ® ^p)Qiy) Q{y){iu, ® QQ{y) , (ii) 

and are thus Lorentz singlets. The appearance of flavor matrices implies that they break 
the SU (4) symmetry, but the contraction of flavor indices between bilinears implies that an 
5*0(4) subgroup survives unbroken. This is the subgroup in which the 4 of 5f/(4) transforms 
in the spinor rep of 50(4): 

Q — ^ (1 ® Ai/2)g , Ai/2 = exp(cj^^eM^) , (12) 

with uj^j_^ an antisymmetric matrix of real parameters. The appearance of operators contain- 
ing bilinears of all possible flavors in eq. (H) implies that no other continuous flavor symmetry 
remains. One can also show that no axial symmetries survive aside from the U{1)a described 
above. 

There is, however, an additional discrete flavor symmetry. All the four-fermion opera- 
tors (of both type A and B) are invariant under the "flavor reflections" generated by the 
transformation^ 

Q^(i®UQ. Q^Q(i®U- (13) 

These transformations may change the sign of a bilinear, 

Q{is ® iF)Q - Q{is ® C,F,)Q = i-f'Qhs ® ^f)Q (14) 

(see App. 0for notation), but this sign cancels in the four-fermion operators since all bilinears 
appear in flavor diagonal pairs. The generators in (p!3D anticommute and square to +1, and 
so generate the Clifford group They lie within the flavor group 5[/(4), but not within 
the 50(4) subgroup (which is composed of matrices containing an even number of ^^). They 
transform as a vector under the flavor 50(4), and thus are combined with this group in a 
semi-direct product as shown in the Table. 

The second stage of symmetry breaking is due to the six operators in Sq^^^\ in which 
the flavor and Lorentz indices are correlated. Two examples appear in App. ^ [eqs. ( p^ ) and 
(|65|)]; here we give another (this time in continuum notation) 

IJ,<U 

- Q{y){i^.u5®i^.^)Q{y) Qiy){i5u^<^^5p)Q{y) ■ (15) 

3ln Ref. we used a different basis for the generators, namely («^5^)- Although the two choices of basis 
are equivalent, that we use here allows a more direct connection to the lattice symmetries. 
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These operators are not invariant under either Lorentz or flavor SO {4) rotations.^] If, how- 
ever, we do a simultaneous Lorentz and flavor hypercubic rotation, then the correlation 
between the indices is preserved. The restriction to hypercubic, rather than continuous. 



rotations occurs because indices are repeated four times [e.g. in eq. (|T5|) ]. The explicit 
form of this diagonal hypercubic rotation is 

Qix) -> (Af/'^^ ® [Af/'^]*)g([A^^]-M , (16) 

with Ay^ an element of the spinor representation restricted to the hypercubic group, and 
A^^ the corresponding member of the vector representation. Since these transformations 
involve simultaneous Lorentz and flavor rotations we refer to the group they form as SVKi^diag- 
As noted above, Sq ^ Ms also invariant under the discrete flavor group the generators 

FF( 

of which transform as a vector under 5'W4 diag- Thus the flavor- rotation group of is 

The symmetry group of 5*6 is given by the intersection of the symmetry groups of its 
components. Combining the above results, we find that r4 x SWi^diag is the symmetry 
group in general, and that in the chiral limit we get the additional U{1)a- Adding back in 
fermion number, charge conjugation and parity, the resulting groups are indeed those of the 
underlying lattice theory [0. This is most easily seen using the presentation of the lattice 
group given in Ref . |jl6| : the r4 is generated by single-site translations, while the SW4^ciiag is 



generated by lattice rotations. 

In summary, we have learned that the bulk of flavor symmetry breaking is caused by the 
subclass of four-fermion operators in which flavor and Lorentz indices are correlated. We 
will now show that the contributions of these operators are suppressed in low momentum 
pion interactions. 



3 Effective Chiral Lagrangian 

In this section we work out the modifications to the chiral Lagrangian implied by the ad- 
dition of the action o^Sq. We begin by recalling the form of the Lagrangian describing the 
long-distance dynamics of the QCD action, 5*4, close to the chiral limit. The spontaneous 
breakdown of SU{4:)i x SU{4:)r to vector SU{4) gives rise to 15 Goldstone modes, described 
by fields 0j. These can be collected into an SU{4) matrix 

15 

E(x) = exp(z0//) , = ^0„Ta, (17) 

a=l 

where / is the pion decay constant (normalized to = 132 Me V). We adopt an unconven- 
tional normalization for the generators, 

Ta = {^M' ^^m5' ^^m^' ^5} • (18) 
Under SUl{4:) x SUr^A), S transforms as 

S(x) — > L^{x)R^ (19) 

*The extra subtractions in the last two compound operators in eq. (|^) are included to remove a Lorentz 
singlet component. 



6 



where L{x) G SUl{4:) and R{x) G SUr{4:). The chiral Lagrangian, correct up to quadratic 
order in meson masses and momenta, is 

4 = y Tr(a^S9^St) - 1/i m fTr{J: + St) , (20) 

with /i a constant of O(Aqcd)- Expanding this in powers of one finds 15 degenerate pions 
with masses given by 

ml = 2/xm [1 + 0(m/AQCD)] • (21) 

The leading order term is the tree-level result, while the corrections come from loop diagrams 
and from higher order terms in the chiral Lagrangian. 

The addition of o^Sq breaks chiral symmetry and lifts the degeneracy of the pions. Gener- 
ically, eq. (|2T|) becomes 



= ca 



A^CD + 2/im 1 + 0(m/ Aqcd) + 0{a' A'^^d) (22) 



where c is a constant of order unity. Contributions proportional to are due to 5*6, and lead 
to massive pions even in the chiral limit. The only exception is the pion with flavor ^5 which 
remains massless because Sq respects the U{1)a symmetry when m ^ 0. It is instructive to 
rewrite the general form as 

= a^A^cD + T + «'A^cDT + + • • • , (23) 

^^QCD ^iQCD Aqcd ^qcd 

where we have dropped all constants. This shows explicitly that we are doing a joint ex- 
pansion in and m about the continuum and chiral limits. If we treat both parameters 
as small, then we can ignore all but the first two terms. The key observation is then the 
following. While the combined effect of all the terms proportional to (and any power of m) 
in eq. ( ^5] ) is to break the chiral and Lorentz symmetries down to the lattice subgroup (since 
5*6 itself breaks the symmetry in this way), it is possible that the term which remains in 
the chiral limit has a larger symmetry. This turns out to be the case, and implies certain 
exact degeneracies between pions in different lattice irreps in the chiral limit. 

To demonstrate this we need to determine the mapping of the operators in 5*6 into the 
chiral Lagrangian. This is done by matching the transformation properties of operators 
under Lorentz and chiral symmetries. We denote the resulting contributions to the total 
chiral Lagrangian i.e. 

£^ = 4 + a24 + 0(a^). (24) 

A precise statement of our result is that the part of £^ without derivatives, i.e. the potential 
V^, is invariant under flavor SO (4) transformations.^] Terms in £^ involving derivatives do 
break 50(4) down to the lattice symmetry group. But for very low momentum interactions, 
such that ~ ml, these terms give contributions suppressed relative to those of by 
m^/AqQj-,. For example in the expression for m^, eq. (|23|) , gives rise to the term 
on the r.h.s, while the terms in £^ containing derivatives give rise to the non-leading term 

X. 

proportional to a^m and a^. Note, however, that for pion interactions with ~ ^qcd ^^e 
5*0(4) symmetry will be broken. 

In the following we work through the operators in Sq matching them into with the 
aim of finding those which contribute to the potential V^. 

A. 

^We adopt the sign convention that — V® plus terms containing derivatives. 
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3.1 Matching of operators in ^f^"^ and Sq^^^^ 

As noted above, most of these operators have the same symmetries as 5*4, and in particular 
do not break the SU{4:) flavor symmetry. They lead to a renormalization of the terms in 
the continuum chiral Lagrangian by a factor of the form 1 + a^Aqf-.^, + a?m^ . They do 
contribute to V^, but their leading contribution to is of order 0{a'^m), and thus of second 
order in our joint expansion. 

There are also terms in S'f''^'^ and Sg''" which are not invariant under rotations. These 
map into operators which themselves are not rotationally invariant, which requires at least 
four derivatives, e.g. 

ETr(aJsaJSt). (25) 

These lead to rotation non-invariance in the pion propagator, and thus direction dependence 
of the extracted pion mass, but they do not break the flavor symmetry. Their contributions 
to are of size a^m^ and thus proportional to a^m?, a^m or a^, and so are of third order 
in our joint expansion. 

3.2 Matching of operators in S^^^^"^ 

The flrst flavor-breaking contribution to the potential comes from the operators in S^^'^^\ 
These can match onto operators without derivatives because they are singlets under rota- 
tions. The explicit mapping is worked out in App. At this stage, the potential retains an 
SO {A) flavor symmetry, since this is the symmetry of all the operators in S^^^^\ 
Using the results of App. P, we flnd that the contribution to the potential is 

+ C2i[Tr(S2)-Tr(e5Ee5E) + /i.c.' . 

+ C4|E[Tr(USe5.S) + /i.c.] 

V 
V 

-vC^Y^Mi^v^ivyX^^)- (26) 

iJ,<V 

The six unknown coefficients Ci are all of size Aq(-;j3. All terms respect the enlarged flavor 
symmetry and are Lorentz singlets. To obtain the result (pGf) , we have collected the contri- 
butions from the different operators in eqs. ( |85| - p7| ), (pl|), (p^), (pG]), and (0), and then used 
Fierz transformations to convert all operators to a form containing a single flavor trace. 

3.3 Matching of operators in ^g^*^^^ 

Our flnal task is to match the operators in i.e. those in which spin and flavor are 

correlated. These are the operators which break flavor down to the lattice subgroup. The 
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key point in this matching is that all the operators transform non-trivially under rotations. 
To represent these operators in the chiral Lagrangian requires the introduction of derivatives, 
since these are the only objects available to introduce non-trivial rotation dependence. This 
implies that the operators in ^^^^ make no contribution to the potential. Thus the result 
eq. (|26| ) is the full result for V^, and is 5*0(4) flavor symmetric as claimed. This is the central 
result of this paper. 

Although complete flavor breaking does not occur at first order in our joint expansion, 
it is interesting to see how it enters at next order. An example of the operators that appear 
when matching S'g ^""^^ is 

[T,y<V,]—^ Y: Tr((9^-9.)Sep)Tr((9,-9.)Sep) . (27) 

This leads to a direction dependent contribution to pion masses that is correlated with their 
flavors. For example, this term renormalizes the mass-squared of a p = pion with flavor C,j 
(j = 1 — 3) by a factor of 1 + a'^AQ^j^, while it does not effect the mass of the pion with flavor 
^4. Since in leading order in our joint expansion, the pion mass-squared is proportional to 
and m, this example shows that the second order correction terms proportional to both 
and a^m in eq. (|23|) break the symmetry down to the lattice subgroup. 



3.4 Contributions of O(a^) 

We close this section with an observation concerning higher order discretization errors. The 
numerical results discussed in the following section indicate that the SO {4) symmetry is 
only weakly broken even at fairly large lattice spacing. This raises the question of whether 
higher order discretization errors break the flavor SO {A). In particular, is the enlarged flavor 
symmetry respected by higher order discretization errors in the chiral limit? 

The analysis of the previous subsection has, in fact, answered this question in the neg- 
ative. Terms in the chiral Lagrangian involving derivatives do break 5*0(4), and, through 
wave-function renormalization, give flavor-breaking contributions of O(a^) to m^. Neverthe- 
less, it is interesting to know whether the "direct" contributions of 0{a'^), i.e. those in the 
potentials V^^^ break 50(4). 

To study this question, all we need to know is that, in the chiral limit, any operator in 
yn>8 jjiust be symmetric under 

f/(l)AX(r4X5iy4,diag), (28) 

for this is the symmetry of the lattice theory for m = 0. The issue is whether all operators 
in V^-^ that respect this symmetry also respect a larger flavor symmetry (as was the case 
for V^). That the answer is negative is shown by the example 

^Tr(e.Se.Se.Se.S) + /i.c. (29) 

u 

This does respect a larger symmetry than (pSj), because (when integrated over all space-time) 
it is invariant under rotations. But its flavor symmetry is the hypercubic group SW4,, and 
not 50(4). Thus flavor is completely broken by higher order terms in the potential. 
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Nevertheless, it is amusing to note that the contributions of such terms to the pion masses 
do respect SO {4). This is because, when calculating pion masses, two of the four factors of S 
are replaced by the identity, and the operator "contracts" to one which is 50(4) symmetric. 

This analysis does not tell us whether the new term (|29|) is proportional to or a higher 
power of a^. To determine this one must study the quark operators of dimension 8 allowed 
by the lattice symmetries, a task we have not undertaken. 

4 Results for pion masses 

The partial restoration of continuum symmetries in applies to all the multi-pion interac- 
tions encoded in the chiral Lagrangian. We consider here only the simplest application of 
this result, namely to the masses of the pions. 

In the continuum, the pions form a 15-plet of flavor SU{A), and are degenerate. On 
the lattice, states are classified by the symmetries of the transfer matrix (a subgroup of the 
symmetries of the lattice action). As shown by Golterman [0, the pions fall into 7 irreps of 
this group: four 3-dimensional reps with flavors ^i5, ^ij and ^j4, and three 1-dimensional 
reps with flavors ^4, ^45 and ^5. Here we have chosen the 4'th direction as Euclidean time. 
These results hold irrespective of the quark mass. 

Close to both the chiral and continuum limits, the pion masses are given by 

M^{Taf = 2fim + a^A{Ta) + O(a'm) + O(a^) , (30) 

with A(Ta) ~ Aqcj-, arising from V^. Since respects flavor 5*0(4) (as well as a discrete 
group which plays no role in relating pion masses), the pions fall into 50(4) representations: 

• A one-dimensional irrep with flavor ,^5, which receives no mass from any of the terms 
in C^, because of the exact axial U{1) symmetry: 

A(e5) = 0. (31) 

• Two four-dimensional irreps 

A(e.) = 
A(U) = 

• One six-dimensional irrep with flavor C,fiu, fJ' ^ ^, with mass 

A(e^.) = ^(203 + 2O4 + 406) (34) 

The mass shifts for these representations are independent, and thus we can make no predic- 
tions for the ordering or splittings. We also cannot predict the sign of the shifts, although 
numerical data indicates that they are all positive with present forms of the staggered action. 



with flavors and ^^5. These receive masses 
— (d + 02 + 03 + 304 + 05 + 306) , (32) 
^(Oi + O2 + 3O3 + O4 - O5 + 306) . (33) 
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Figure 1: Splittings between non-Goldstone and Goldstone pions at /? = 6 in the quenched 
approximation Results are in lattice units. The solid horizontal lines show the average 
value for the three pairs which are predicted to become degenerate in the chiral limit. Errors 
are approximate, as they have been obtained ignoring the error in the Goldstone pion mass, 
and the correlation between the masses of pions in different representations. 

We can finally see the degeneracies predicted in the chiral limit at finite lattice spacing. 
The lattice irreps with flavors and ^4 join together, as do those with flavors and ^45, 
and those with flavors and ^44. 

How well do these predictions compare to numerical data? Two groups have calculated 
the masses of all seven pion irreps with errors small enough to allow a test. Ishizuka et al. 
made an extensive study of flavor symmetry breaking several years ago . They worked in 
the quenched approximation at /? = 6.0, and used two quark masses m = 0.01 and m = 0.02. 
These masses corresponding roughly to m^/S and 2m<(/3, with rris the physical strange 
quark mass, and thus are moderately small. The chiral expansion parameter (m^/47r/^)^ 
is about 25% for such masses. The second calculation has been done recently by Orginos 
and Toussaint |Q. They have done extensive calculations using a variety of perturbatively 
improved staggered fermion actions using moderately light dynamical quarks. 

It turns out that the predicted degeneracies are observed in almost all the data sets, even 
at the largest quark masses. This is illustrated in Figs. and |^, which show, respectively, the 
results of Ref. and a subset of the results from Ref. [Q. For each of the six non-Goldstone 
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Figure 2: Splittings between non-Goldstone and Goldstone pions in dynamical simula- 
tions 1^. Notation is as in Fig. |^. The upper two sets of points come from simulations 
with unimproved gauge and fermion actions, while the lower two use improved gauge and 
fermion actions (the latter being the "fat+Naik" action). For clarity, some points have been 
offset horizontally. 



pion irreps, we plot the difference 



M^iTaY - M^{^,y =a^A{Ta) 1 + 0(m/AQCD) + 0(a^A4cD 



(35) 



which removes the common contribution proportional to the quark mass. We order the flavors 
so that pairs which are predicted to become degenerate in the chiral limit are adjacent. For 
all pairs except one (^j and ^4 at the heaviest quark mass with unimproved fermions in Fig. |^), 
degeneracy is observed within errors. n Furthermore, there is, in most statistically 
significant difference between the masses of the pions in different 5*0(4) representations. 

The numerical data are thus consistent with a flavor SO (4) symmetry, not only in the 
chiral limit but also for non-zero quark masses. As noted in the previous section, this 
symmetry is broken by both the higher order terms in eq. (|35|). In particular, SO (4) breaking 
by the correction proportional to m would be twice as large for the heavier masses in each 

^In simulations, the masses of the various pions are correlated, and thus the differences between them are 
more significant than the error bars suggest. It would be interesting to directly calculate the errors on these 
differences. 
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of the data sets. It is apparent from the figures that such S0{4) breaking is small, and thus 
that the bulk of this correction is 5*0(4) symmetric. This result is not explained by the 
chiral Lagrangian analysis. 

It is interesting to compare the two data sets in Fig. ^. The two sets are approximately 
matched in the following sense: for both the smaller and the larger quark masses, the values 
of aMp and aM^(^5) from the unimproved and improved simulations agree closely. This 
means that they have roughly the same lattice spacing and quark mass, and can thus be 
used to observe the effect of improvement on flavor symmetry breaking. At non-zero quark 
mass the flavor breaking is substantially reduced, particularly for the pions with flavor ^^5. 
On the other hand, if one extrapolates linearly to the chiral limit, the effect of improvement 
is much smaller. It is also noteworthy that, for improved fermions, for which one can do 
the chiral extrapolation with little uncertainty, the pattern of flavor breaking is consistent 
with only the C4 term in eqs. ( p^ - |3^) being substantial. We have no explanation for this 
observation. 

Finally, we can use the data to look at the absolute size of the discretization errors in 
the chiral limit. The quantity shown in the Figures should be of size (aA)^, with A ~ Aqcd- 
Performing a linear extrapolation to m = 0, and taking the result for the flavors as 
representative, we find the correction to be 0.031 for the quenched data, and 0.096 (resp. 
0.11) for the improved (resp. unimproved) dynamical simulations. The lattice spacing for 
these simulations (obtained from extrapolating the rho mass to the chiral limit) is roughly 
1/a = 2 GeV for the quenched data, and 1.3 GeV for both dynamical simulations. Thus 
one finds that the scale describing the discretization errors is A ^ 0.8, 0.7 and 0.7 GeV for 
the three cases. These rather large values illustrate the need for further improvement of 
staggered fermions. 



5 Non-Goldstone particles 

An interesting question is whether the partial symmetry restoration found in the pion sector 
extends to other particles such as the vector mesons and light baryons. A naive argument 
against such an extension goes as follows. Even for pions, the enlarged symmetry is not 
respected by terms proportional to a^p'^. For other hadrons, e.g. the p meson, one has 
Ip^I = rn?p 3> m^, and so the flavor-breaking terms are not suppressed compared to other 
discretization errors. We have studied this question in detail for vector mesons, and find, as 
explained in this section, that indeed symmetry restoration does not occur. 

The essential difference from pions is that the vectors remain massive in the combined 
chiral and continuum limits. Indeed, the simplest way to include them in the chiral La- 
grangian is by expanding in inverse powers of their masses [0. At leading order they are 



sources having fixed velocity. Higher order terms fit naturally into a chiral expansion because 
the expansion parameters match: l/nip ^ 1 / {Att f.,^) . 

The key point is that there is no obstacle to constructing terms, at leading order in the 
chiral expansion, with the same transformation properties as those contained in Cq^^^\ In 
particular, unlike for pions, one does not have to build Lorentz non-singlets using derivatives. 
Instead, one has two new Lorentz vectors, namely the rho field itself pp and the four-velocity 
of the heavy source, Vp. What is particularly important is that additional factors of Vp can be 
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added at no cost in the chiral expansion. To see how this works in detail is straightforward but 
tedious. We give only an example. The result of matching the operators ([V x l^] — [A x A]) 
and ([V^ X V^] — [A^ x A^]) onto operators in the chiral Lagrangian includes the following 
mass terms, 

MpUupf,^ + J2 MpUuPu^ + Mp^upf^^v^v^ , (36) 

with unknown coefficients. Here is a flavor matrix containing the SU{A) 15-plet of vector 
mesons. There are also similar terms with "axial" flavor. Note that both the second and 
third terms break Lorentz invariance and flavor symmetry. 

Performing a similar analysis with all the operators in we find the following form for 
the a? contribution to the vector meson mass matrix in the chiral limit: 

ci = RiY.Mpip,] + R2T.Mp%pM 

fj. fj. 

+R3 Mpl^uPt,^ +RiYl MpUu5Pf.^u5] +Rrj Mp%xp^^ux] 

+Re Mp%Pf^^ +R7Y Tr[pl^M5pM^M5] + RsY Mp%<^Pf^^,^'^] 

+Rg Y MpU^Pf^^^ + ^10 Y Tr[pi^45P/.^45] + ^11 Y Mpi^iuP^,Uu] ■ (37) 

Here we have specialized to the rest frame, in which v = and = 0. The unknown 
coefficients Ri are of size o^Aqq^)- These eleven coefficients break the vector meson octet 
down into the eleven lattice irreps found by Golterman |T^. These are the states created by 
the bilinears 

Q(7fc®C5)Q, Q(7fc®^m)Q, Q{ik®U)Q, 

Q{lk®ik)Q, Q{lk® irab)Q, Q{lk® iib)Q, 
Q{lk ® ikb)Q, Q{lk ® iraA)Q, Q{lk ® ilm)Q, 

Q{-fk ® ^k4)Q, Qilk ® ^mk)Q, (38) 

where k ^ I ^ m. Thus there are no predicted degeneracies. Note that the terms in the 
last line of eq. (|37D , which are present because of the vector f^, are necessary to lift all 
degeneracies. 

The conclusion we draw from this example is that the pions are a special case, and that for 
all other particles there is no reason to expect a partial restoration of the lattice symmetry. 
At the present time there are no numerical results accurate or extensive enough to test this 
prediction. 



6 An "Aoki-phase" for staggered fermions? 

As the quark mass is lowered, there comes a point at which the effects of the quark mass 
term and the leading discretization errors are comparable. This occurs for physical masses 
m ~ a^, implying lattice masses am ~ a^. At this point the competition between the two 
contributions can lead to unexpected patterns of symmetry breaking. For Wilson fermions 
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one has the possibihty of spontaneous breakdown of flavor and parity symmetries, i.e. an 
"Aoki-phase" containing exactly massless Goldstone bosons 0. What we wish to point out 
here is that an analogous phenomenon can occur with staggered fermions — the 5*0(4) flavor 
symmetry of the terms can be broken spontaneously. The method we use is a direct 
generalization of the work of Ref. |]^ on the Aoki-phase for Wilson fermions. 



That symmetry breaking can occur is clear from the form for pion masses, eq. (|30|). If the 
corrections of O(a^) (called A(T(j) in the equation) are negative for some pion flavors, then 
MT^iTaY will become negative for small enough quark mass. This signals vacuum instability 
and the possibility of spontaneous symmetry breaking. As noted in Sec. H, present versions 
of the staggered action appear to have all A(Ta) positive, in which case symmetry breaking 
does not occur.0 Nevertheless, it is possible that alternate improved versions will lead to 
symmetry breaking, and we discuss this possibility briefly. 

Before doing so, however, we stress an important difference between the Aoki-phase 
with Wilson fermions and the present analysis. With Wilson fermions, the flavor symmetry 
in question is an exact symmetry of the lattice theory, and its breaking leads to massless 
Goldstone bosons at non-zero lattice spacing. In the present case, the flavor symmetry 
is only approximate — it is violated by terms of size a^. These terms lead only to small 
corrections to the masses of most of the pions in the region of interest (M^ ~ a^A^). The 
only exceptions are the Goldstone pions (0g below), which become pseudo-Goldstone bosons 
with Ml ~ a^A^ 



To simplify the discussion we assume that only the C4 term in £^ is present [see eq. (p6|)] 



and assume that C4 itself is negative. A straightforward calculation flnds the following 
features as the magnitude of the quark mass is reduced. For |m| greater than the critical 
value iTLc = 24a^|C4|/(/i/^) ~ a^A^, the condensate is given by (E) = sign(m)l, and the mass 
formulae of Sec. ^ apply with m — |m|. When |m| = rric the pions with flavor become 
massless [see eqs .(|31|)-(|3^)], and there is a second order phase transition. For |m| < it 
is useful to introduce the angle 6 deflned by cos^^ = m/mc, with < < vr. Minimizing the 
potential, one flnds that the condensate swings from +1 to —1 according to 

(S) = exp(i6'[^^ ■ n^]) = cos6' + isin6'[^^ ■ n^] , (39) 

where is an arbitrary vector of unit length. For < < tt, flavor 5*0(4) is broken to an 
5*0(3) subgroup, and there are three Goldstone pions. As noted above, these are not exactly 
massless because terms in of 0{o?m) and O(a^) break the 50(4) explicitly. 

It is instructive to work out the full spectrum of pions in the broken phase. We display 
these for = 5^4, i.e. for the case in which the condensate is (E) = exp(i6',^4). The 
pions are labeled by their flavor, deflned by the direction of excitation of the condensate, 
E = (E) exp(z(/)aTa). The results are given in units of 16a^|C*4|//^. 

• Flavor Ta = i^^j, with j = 1, 2, 3. For each j, these two flavors mix, and there are 
two mass eigenstates: 

<Pg = cos ^0(e,) + sin ecl){C4j) , = ; (40) 

(Pngi = cos ^0(^4^) - sin 9(j){Q , M^^i = 1 • (41) 



^When m becomes negative, the vacuum expectation value (S) flips from +1 to —1, and the result ( pOj ) 
remains valid with m — > Iml. 
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At the end-points, the Goldstone modes are pure flavor ^j, but at the mid-point they 
become pure flavor ^4^. 

• Flavor Ta — ^i- This single state has 

M^^Ssin^^. (42) 

This vanishes at the end-points, due to the restoration of SO {A) symmetry, which 
makes the masses of all the flavor states equal. 

• Flavor Ta = i^jk- Three states with 

= cos^ e . (43) 

These become degenerate with the other tensors at the end-points, due to 5'0(4) 
restoration. This degeneracy recurs at the mid-point, where all six tensors are massless. 
This is not, however, caused by 5*0(4) restoration: the .^4^ states are massless due to 
the spontaneous breaking of flavor 5*0(4), while the C,jk states are massless because of 
an accidental axial 5*0(3) symmetry of the potential about its classical minimum at 
m — 0. The latter states will presumably become massive when quantum corrections 
are included. 

• Flavors = ^{^5. Three states which are always massive: 

= 1 + cos^ e . (44) 

These become degenerate with the at the mid-point, a result that can be understood 
as due to the restoration of the U{1)a symmetry. 

• Flavors = (^5,^^45. These mix according to 

(f>NG2 = cos ^0(^45) + sin ^0(^5) , M^a^ = 2 ; (45) 

(f)NG3 = cos ^0(^5) - sin ^0(^45) , M^G3 = 3 . (46) 

The degeneracy between flavors ^45 and ^4 at the mid-point is again due to U{1)a 
restoration. 

7 Conclusions 

We have studied discretization errors for staggered fermions, by constructing a sequence of 
effective Lagrangians. Our central result is that, close to the continuum limit, the potential 
term in the chiral Lagrangian describing pion interactions respects a chiral symmetry group, 
U{1)a X [r4 XI 5*0(4)], which is much larger than that of the underlying lattice theory. Full 
5*0(4) rotation invariance is also maintained. The enlarged symmetry predicts degeneracies 
in the pion masses which are observed in numerical data with surprisingly small deviations. 

While this result is pleasing, it has limited applicability. The enlarged symmetry is 
broken in pion interactions by terms of O(a^), O(a^m) and O(a^p^), and also does not hold 
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for hadrons other than pions. We do not know whether it holds for other pion properties, 
e.g. the decay constants. 

An interesting question is whether our result can be used to simplify the improvement 
of the staggered fermion action. Complete non-perturbative improvement at O(a^) is pro- 
hibitively difficult because of the number of additional operators that are needed. If flavor 
symmetry were partially restored, if only approximately as at finite quark mass, then one 
should be able to get away with fewer operators. An extreme example is the fiavor breaking in 
pion masses alone. After symmetry restoration there are only three independent differences, 
so one could imagine a non-perturbative tuning of the coefficients of just three operators so 
as to eliminate these differences. This approach was suggested by Lepage 0, who pointed 
out that tree level improvement of the staggered action contains only three-flavor breaking 
terms (which can be written either as four-fermion operators or as bilinears with fat- 
links 0]). One can show that these three operators do give independent contributions to the 
pion splittings, so that such tuning is feasible. On the other hand, the flavor-breaking part 
of the potential contains six independent parameters (Ci — Cq in our notation), and so to 
remove it non-perturbatively would require introducing six improvement operators. Further- 
more, there is no simplification if one wishes to completely improve the full spectrum, for 
which there is no symmetry restoration. Thus it is unclear to what extent the limited fiavor 
restoration will help. Numerical tests are needed to investigate this issue, and exploratory 
work has been done in Refs. P, ||, ^ . 
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A Enumerating operators of dimension six 

In this appendix we explain how to determine all operators of dimension six which are singlets 
under the lattice symmetry group. These are the operators which must be added to remove 
O(a^) errors from the staggered action (modulo possible redundancies). Conversely, when 
converted to continuum form, they are the operators which appear in the continuum effective 
Lagrangian at O(a^). Luo has cataloged these operators previously 0, but did not find all 
of them. 



We use the notation of Refs. |T^, ^ |2T[], the relevant parts of which we briefiy review 
here. The continuum theory describes four fermion fiavors, which are collected into a field 
Qa,a, with a and a respectively spinor and fiavor indices. Bilinears are labeled by a spin and 
a fiavor matrix 

Qa,abs)al3{^F)abQl3,b = Qa,ahs ® ^f) aa,f3bQ I3,b = Q{lS ® ^f)Q ■ (47) 
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The spin matrices are labeled by a hypercube vector S^, with components or 1, 

7s = 7f 7f^7f^7f • (48) 

The flavor matrices are labeled similarly by the hypercube vector F^, except that they are 
built out of the complex conjugate matrices, = 7*. This is simply a convention, since 
these two bases are unitarily equivalent. Finally, we use abbreviations such as 

7^i>y = li^lu , 7/^5 = 7m75 • (49) 



The lattice bilinears which correspond in the continuum limit to (^TQ are 



J2xiy)ciis®^F)cDxiy)D, (50) 

CD 

where y labels the 2* hypercubes, C and D are hypercube vectors, and the hypercube field 
is defined in terms of the underlying staggered fermion by 

x{y)c = lx{y + C). (51) 



The matrices (75 ® ^p) are unitarily equivalent to (75 ® ^f)- 

It will be useful in the following to define hypercube fields at zero physical momentum, 
which are obtained by averaging over all A^^ hypercubes: 

xc = ^Y.xiy)c. (52) 

The 16 xc transform in the defining representation of the lattice symmetry group. We always 
use the argument y when referring to the fields which live on individual hypercubes, and 
drop this argument for the zero-momentum fields. 

We begin with a comment on fermion bilinears of dimension six. Luo lists seven such 
operators. We note however that one linear combination, which in Luo's notation is 

x{y){V'p-V>V')x{y) (53) 

has negative parity under lattice charge conjugation and is thus forbidden. Furthermore, 
there are two additional operators, not included explicitly by Luo, 

m'^xiyWxiv) and m^xiy)xiy) ■ (54) 

These are, however, redundant, since they can be absorbed by changing the normalization 
of the fields and quark mass. 

We now turn to four-fermion operators. We want to construct all such lattice operators 
which are singlets under the lattice symmetry group and which correspond to operators 
of dimension six in the continuum. The latter requirement excludes operators containing 
derivatives or factors of the quark mass. Since the quark mass does not appear, the con- 
struction is effectively in the chiral limit, and thus the lattice operators must also be singlets 
under U{1)a- 

We construct operators by the following steps. These make use of the fact that, to 
generate the lattice symmetry group, one can replace rotations about a point on the lattice 
with those about the center of a hypercube, and similarly for reflections. These two choices 
differ only by translations. 
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1. Multiply two bilinears residing on the same hypercube, and sum over hypercubes. 



o{s, F, s\ F') = J2xiy)iis ® ^F)xiy) xivKis' ® ^fMv) ■ (55) 

y 

This operator is manifestly invariant under translations by two units. Gauge invariance 
is maintained (in a way which maintains rotation and reflection properties) by including 
the usual average of products of gauge links along the shortest paths between the x 
and X- This can be done in two ways: joining the x X within the bilinears, or 
between the bilinears. If we ignore the gauge links, these two choices are related by a 
Fierz transformation. This remains true in the presence of the links, up to additional 
higher dimensional operators involving gauge fields. Thus we need consider only one 
linear combination of color structures. We do not need to specify our choice for the 
following analysis. Luo uses the linear combination which arises in perturbation theory 

I (contract within bilinears) — | (contract between bilinears) . (56) 

When we compare our operators to his, we are implicitly using this choice. 

2. Form linear combinations of the 0{S, F, S', F') which are singlets under the symmetry 
group of the hypercube — rotations and reflections — and also charge conjugation. We 
label this group . The transformation properties of bilinears under have been 
worked out by Verstegen and are listed in his Tables 3 and 4. The bilinears fall into 



56 irreps. To form a four-fermion operator which is a singlet one must combine two 
bilinears which live in the same irrep. Since most irreps appear multiple times, 
there large number of ways of doing this. 

3. At this stage the operators are not invariant under single site translations. This can 
be accomplished by applying the projector 

ni(i+^.)' (57) 

with 7^ the translation operator in the /x'th direction. Note that the 7^ commute when 
acting on operators with even fermion number, such as the four-fermion operators under 
consideration. Thus the order of the factors in the product in (|57|) is unimportant. 



4. The resulting operators are singlets under all lattice symmetries except U{1)a- We 
now apply the constraint that they do not contain derivatives. This can be done by 
taking the tree level matrix elements with all external quark states having zero physical 
momentum, i.e. mod(p^a, vr) = 0. This matrix element will vanish if the lattice 
operator corresponds to a continuum operator containing one or more derivatives. But 
in evaluating this matrix element each of the hypercube fields can be replaced by 
its zero- momentum counterparts: x(2/)a Xa- Thus we are led to consider linear 
combinations of operators of the form 



0'{S, S\ F, F') = n + Vxils ® ^F)xxils' ® ^F')X (5^ 
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We can simplify this using the known translation properties of the zero-momentum 
hypercube field 



%X={I®^^.)X ^MX(75®eF)x= (-l)^"x(75®eF)x, (59) 

where 

F^ = mod(^F,,2). (60) 

We thus find that 

0'{S,S',F,F') = ni(l + (-)(^-^')-)x(75 ® ^f)xx{1s' ® ^f')x • (61) 

This operator vanishes unless F = F', which implies F = F' . In other words, if the two 
bilinears have different flavor, then the projection onto a translation singlet produces a 
continuum operator containing derivatives. Thus we conclude that, in step 2, we must 
keep only those bilinears in which F = F' . 

5. Finally, we select those from the resulting list which are singlets under U{1)a- 

Our claim is that this procedure produces all singlet dimension-6 operators. 

Thus we must determine which of the lattice four-fermion operators satisfying F = F' 
("diagonal in flavor") are singlets under hypercube rotations and reflections. This is a 
straightforward, though tedious, group theoretical exercise which can be done using the 
tools presented in p2[. First, one finds that to make a singlet it is necessary that 5* = 5*', 
so that the operator is diagonal in both spin and flavor. Then the problem reduces to 
picking out the singlets in the product of vectors in which each index appears twice, e.g. 
7/i ® 7^4 ® ® Ct^- Note that the "square" of a flavor matrix transforms in the same way as 
the "square" of a spin matrix. The result is that there are 35 singlets |pi|| : 

A) There are 25 operators in which the spin and flavor indices are separately contracted. 
For the spin matrices we use the notation 

S=l^l, P = 75®75, V^ = ^7^®7^, 

^ = 7/^5 ® 75m , T = J2lf^'^^ lyf^ ' (62) 

while for flavor matrices we use the same notation with 7^ —>■ The 25 operators are 
simply the products of the 5 possible spin structures with the 5 flavor structures. The 
notation we use is exemplified by 



X = E E x(y)(7M5 ® ^up)x{y) x(y)(75M ® U)x(y) • (63) 
B) The remaining 10 operators have the spin and flavor matrices coupled. They are 



K X Vf^] = T.x{y){i, ® ^My) x{y)hp ® Ux{y) , (64) 



20 



along with [V^ x A^], [A^j x V^] and [A^ x A^] defined analogously; 



K^T^] = J2xiy)il^^'S)^My)x{y)il^^^U)xiy) 



[S X A] 


= ^9, [SxV] = 


^ 13 ! 






[A X S] 


= -^3, [yxS] 


= -^11 ) 






[P X V] 


= ^10, [PxA] = 


= 17 5 






[V X P] 


= ^4, [AXP] = 


—-^ 15 5 






[TxV] 


= ~(|-^ 3 + -^le) 7 


[T X A] = 




+ 12) , 


[VxT] 


= — (1-^7 + -^is) 5 


[A X T] = 


(-^8 


+ 14) , 


[T, X V,] 


= (~|-^ 3 + -^le) 5 


X A,] 


= (- 


+ \ 12 


K X T,] 


= ("l-^ 7 + -^is) ) 


[\ X T^] 


= (- 


'•^8 + \^ 14 



- (7/. ® C/.i^5)x(l/) (7m ® i^vi,)x{y) , (65) 

with [T^ X V^], [y4^ X Ty\ and [T^ x A^] defined analogously; and finally [T+ x T+] 
and [T_ x T_], the definitions of which we do not reproduce since they are ?7(1)a 
non-singlets. 

Next we select from these operators those linear combinations that are \J{\)a singlets. 
This is automatically true for those operators consisting of odd bilinears.0 There are 16 of 
these, all of which were found by Luo. The relation to Luo's basis is: 

(66) 
(67) 
(68) 
(69) 
(70) 
(71) 
(72) 
(73) 

To determine which of the operators composed of even bilinears are \J(\)a invariant, we 
note that any such operator must, under fiavor-spin Fierz transformation, transform into an 
operator composed of two odd bilinears. Conversely, if we Fierz transform each of the above 
16 operators, we will find all even-even diagonal operators. Note that some of the odd-odd 
operators Fierz transform back into odd-odd operators, so that we can end up with less than 
16 even-even operators. Using the Fierz tables collected in App. A of Ref. [^, we find that 
there are 8 even-even \J{\)a singlets: 

{[5x5]-[PxP]} = i(J-i+^2), (74) 

{\V,xV^,-\A^xA^W = i(J-i-^2), (75) 

{[SxP]-[PxS]} , (76) 

{[V, X A,] - [A, xV,]}, (77) 

{[VxV]-[AxA]} , (78) 

{[VxA]-[AxV]} , (79) 

{[SxT]-[PxT]} , (80) 

{[TxS]-[TxP]} . (81) 

The first two are related to two of Luo's operators, while the remaining six are new. Note 
that no [T x T] operators are U{1)a singlets. 

^An even/odd bilinear is one in which the x s-nd x s-re separated by an even/odd number of hnks. 
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The continuum versions of these operators (which appear in the effective Lagrangian for 
unimproved or partially improved staggered fermions) can be obtained simply by replacing 
the x's with Q's, and removing the bar on the matrices. In the text we use the same notation 
for the corresponding continuum operators, e.g. [5x1^], as we do for the lattice operators. 



B Mapping of operators in Sq^^^^ into 

The operators in S^^^^^ break into three classes according to the spin of the bilinears. 
B.l Operators with spin structure V or A 

The general form of these operators is a linear combination of [V x F] and [A x F] , with F 
denoting the flavor. The flavor structure of these operators is a sum of terms of the form 

Op = ±Y^ (Q^{^^ (g) Fr)Qr ± Ql(7^ ® (82) 

where F^'^ are hermitian matrices in flavor space. Both ± signs are positive (resp. negative) 
if the spin is V (resp. A). For the axial operators, the internal minus sign is due to the extra 
75, while the external sign appears because for spin A the Dirac matrix is i'y^'j^ so as to be 
hermitian. To determine the form of the corresponding operator in we promote Fl^r to 
spurion fields. If they transform as Fl LFlL^ and Fr — > RFfjR\ then Op is invariant 
under chiral transformations. We then build all the operators out of S which are chiral 
singlets, quadratic in F (linear and cubic terms are forbidden by the F —F symmetry, 
while quartic terms would be of 0(a^)), and parity invariant [F^ Fr and S S"''). At 
the end, we set F^ = Fr = F, where F is the flavor matrix appearing in the operator under 
consideration (e.g. ^5 for [V x P]). 

It turns out that there is only one non-trivial operator into which Op can map 

Of ^ cTr(FiSF^St) . (83) 

The constant, c, is unknown, but has the same magnitude (and sign) for both V and A spins. 
The sign is the same because the overall sign in ([8^ ) cancels with the "internal" sign coming 
from the fact that this operator contains one factor each of Fl and Fr. A second operator 
allowed by the symmetries, 

Tr(Fi)Tr(EFRSt) + Tr{Fn)ME^ Fl^) , (84) 

turns out to be a field independent constant when one sets F^ = Fr = F. This is either 
because Tr{F) = (true for flavor P and T) or because SS''" = 1 (for flavor singlet bilinears). 
Finally, the operators Tr(F|)Tr(EE''') and Tr(F|ES''') are trivial for any F. 

Thus we find that the operators map into the following terms in the potential V^: 

[V X S] and [Ax S] — > 4c , (85) 
[V X P] and [AxP] — > cTr{^5^^5^^) , (86) 
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X T] and [A X T] cJ^M^^^E^^^E^), (87) 

{[V xV]-[Ax A]) and {[A x - [V' x A]) 

^ cj] {Tr(e,Se,St) - Tr(USe5.St)} . (88) 

The fact that the same constant c appears in all these transcriptions is crucial for the last 
of these mappings, eq. (^8l). A check on this result is that the operator is invariant under 
the axial f/(l) only if the relative coefficient of the two terms is —1. Apart from this, the 
fact that the same constant c appears in all transcriptions does not lead to useful relations 
because the coefficients of the underlying quark operators are different and unknown. 



B.2 Operators with spin structure S or P 

The chiral structure of operators of the form [S x F] and [P x F] is 

o'f = (gZ(i ® Fl)Qr ± g^(i ® Ffi)QL)' , 

with the ± sign corresponding to S or P. Here the spurious must transform as Fl 
and Fr RFrLK The non-trivial operators into which O' maps are 



O' 



±CiTr(F^S)Tr(FiSt) 

+ C2 [Tr(F«S)2 + Tr(FiSt)2^ 

+ C3 IMFr-LFrI:) + Tr(FiStFiSt) 



From this we find that 

[S X V] and [P X V] 

[S X A] and [P x A] 



±ciTr(e,S)Tr(^,St) 

+ C2 [Tr(e.S)2 + /i.e." 

+ C3 [Tr(e.Se,S) + /i.e.] 
±CiETr(aS)Tr(e5.St) 

V 

+ C2E[Tr(e.5S)Tr(e5.S) + /i.c.] 

V 

+ C3E[Tr(e.5Se5.S) + /i.c.] . 



For the operators which are linear combinations, we find 



{\S X S] - [P X P]) 



ci [Tr(S)Tr(St)+Tr(e5S)Tr(e5S^) 
+ C2 [Tr(S)2 - Tr(e5S)2 + /i.e." 
+ C3 Tr(S2) - Tr(e5Se5S) + h.c. 



(89) 



(90) 



(91) 



(92) 



(93) 



with ([S* X P] — [P X S*]) mapping into the same three operators except that the C2 and C3 
terms change sign. A check on these results is that the relative plus sign in the C\ term and 
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minus signs in the other terms are those needed to make these operators invariant under 
U{1)a- Finally, for the flavor T only the Ci term survives 



{[S xT]-[Px T]) 2c, Tr(e^,S)Tr(e,^St) . (94) 
B.3 Operators with spin structure T 

There are two such operators, [T x V] and [T x A]. These have a slightly different chiral 
structure to that of the S and P operators, 

Of = Y1 {Q~l{i^.u ® FL)QRf + (0^(7./. ® Fr)QlY . (95) 

In particular, there are no cross terms between LR and RL bilinears. It follows that the 
mapping is as in eq. (|90|) except that there is no ci term, and the other terms have new 
coefficients and C3. From this we find that 

[TxV] +4^ [Tr(e.S)2 + /i.c." 

+ C3E[Tr(e.Se.S) + /i.c.] , (96) 

V 

[TxA] -4E[Tr(e.5S)' + /i.c." 

V 

+ C3E[Tr(ase5.S) + /i.c.] . (97) 

V 
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